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Abstract. The presence of a kHz frequency in LMXBs has been expected from scaling laws, by analogy with the QPO 
phenomenon in HMXB X-ray pulsars. Interpretation of the two kHz frequencies, observed in accreting neutron stars, in terms 
of non-linear resonance in strong-field gravity led to the prediction of twin QPOs in black hole systems, in a definite frequency 
ratio (such as 2/3). The imprint of a subharmonic of the 401 Hz rotation rate in the frequencies of the QPOs detected in the 
accreting millisecond pulsar is at once a signature of non-linear resonance and of coupling between accretion disk modes and 
the neutron star spin. 



INTRODUCTION 

This is a story of the origin of high-frequency quasi- 
periodic oscillations (HF QPOs) observed in the X-ray 
flux of low-mass X-ray binaries (LMXBs). We believe 
that they really are oscillations, resonant oscillations of a 
body of fluid nearly in equilibrium — the accretion disk. 
Accretion disk are so poorly understood that their de- 
tailed numerical and analytic modeling has failed to re- 
veal phenomena approximating the behavior of real HF 
QPOs. Even so, they must exhibit behavior generally ob- 
served in any sufficiently complicated system, e.g., reso- 
nances. 

Physics has principles and techniques allowing a 
meaningful discussion of the behavior of systems whose 
detailed structure is not known. The application of these 
principles has allowed us to successfully predict yj, |2|] 
millisecond variability in LMXBs, at a time when at most 
~ 50 Hz oscillations had been seen, and the appearance 
of pairs of hHz QPOs with rational (as in 3:2) frequency 
ratios in black holes, at a time when only single QPOs 
had been reported. We see no reason to abandon these 
principles now, when the field of QPO modeling has been 
thrown into turmoil with the discovery Q of a subhar- 
monic frequency difference in the accreting millisecond 
pulsar. This is where they take us. High-frequency QPOs 
are an accretion disk phenomenon reflecting fundamen- 
tal properties of gravity. A spinning neutron star may di- 
rectly perturb the accretion disk and influence the fre- 
quencies of its oscillation via a non-linear coupling. 



AN ACCRETION-DISK PHENOMENON 

There is little doubt that QPOs reflect rapid variations 
of flow in accretion disks. The argument is simple. Es- 
sentially the same phenomenon is seen in three types of 
mass-transfering binary systems. In all cases the donor is 
a low-mass star and the bright source is a compact object 
in state of accretion: a white dwarf, a neutron star, or a 
black hole. The observed time-scale of variation depends 
on which of the three, so the donor canot be responsible 
for (and is not expected to exhibit) the rapid variability. 
The white dwarfs have a surface capable of emission and 
an accretion disk around it. The neutron stars likewise, 
and also an occasional jet. Black holes (microquasars) 
have an accretion disk and a jet but no surface. There- 
fore, if the phenomenon has the same origin (see e.g., |4] 
for a comparative phenomenology) it can only occur in 
the single structure that the three types of systems have 
in common: the accretion disk. 



FUNDAMENTAL OR ACCIDENTAL? 

Some fifteen years ago, when quasiperiodic variability 
of luminosity had already been detected in white dwarfs 
(cataclysmic variables), in strongly magnetized accret- 
ing neutron stars (X-ray pulsars) and in weakly (if at all) 
magnetized accreting neutron star, the (tacit) question of 
the day was whether these were separate phenomena re- 
flecting accidental properties of the systems, such as their 



rotation rates, magnetic fields, and the such, or whether 
instead, this was a phenomenon related to the fundamen- 
tal frequency of the system, (2n)VK = yjGM/r 3 0. 
Much confusion arose, because the highest frequency 
then observed in LMXBs was only ~ 50Hz<< Vk- Ig- 
noring this particular QPO, the authors of ref. 1 1 ] opted 
for the Keplerian model of the sub-Hertz QPO reported 
in the X-ray pulsar EXO 2030 + 375 0, and suggested, 
by analogy, that in a neutron star with magnetic field 
< 10 8 G a frequency close to the orbital frequency in the 
marginally stable orbit of general relativity, Vs:(r ms ) « 
2.2kHz (Mq /M), should be observed. Of course, a fre- 
quency of about 2 kHz per inverse solar mass is precisely 
what has been discovered with RXTE in LMXBs, so it is 
reasonable to examine in detail the line of reasoning. 



SCALING WITH RADIUS 

The first assumption 1 1 ] was that the frequency should 
(more or less) exhibit the fundamental scaling of Ke- 
pler's third law, v °< r~ 3 / 2 . This scaling is roughly fol- 
lowed. Indeed, for X-ray pulsars, where the disk ends 
close to the magnetospheric radius r w 10 8 cm, the fre- 
quency is v ~ 0.2 Hz 1 6] (see also Mark Finger's talk 
in this conference Q]), in white dwarfs with radius r ~ 
10 9 cm the frequencies cluster at v ~ 0.04 cm, and for 
neutron stars in LMXBs, with their stellar radius ~ r ms ~ 
10 6 cm 1 8], the QPO frequency is v — 1 kHz. Note that 
in the first two cases the radius depends on the stellar 
properties. For the LMXB neutron stars, the inner disk 
radius was expected to be determined by properties of 
Einstein's gravity HHO]. 



SCALING WITH MASS 

The second assumption was that once the direct influence 
of the stellar surface and magnetic field is negligible, the 
characteristic frequency is determined by gravity, so that 
the characteristic radius determining the QPO properties 
is proportional to GM/c 2 , and hence Vqpo 1/Af, the 
inverse mass of the neutron star (or a black hole). For an 
elaboration of this point see the companion contribution 
by Abramowicz and Kluzniak 1 10]. 

This scaling is satisfied quite well. For ~ 1M© neutron 
stars, Vqpo ~ 1kHz, for a ~ 10M Q black hole Vqpo ~ 
0.1 kHz, and for the 4 • 1O 6 M black hole at the Galactic 
center the recently reported period (of an infrared flare) is 
17 minutes 1 11], i.e., Vqpo ~ lmHz. It seems quite safe 
to assume that the HF QPOs are closely related to orbital 
frequencies in Einstein's gravity. 



TWO VARIABLE CHARACTERISTIC 
FREQUENCIES 

The puzzle of twin kHz QPOs in neutron stars is this 
Q|: the two frequencies have characteristic values, and 
yet they vary in time (within minutes). It is easy enough 
to find two characteristic frequencies for neutron stars 
or black holes, e.g., the maximum radial epicyclic fre- 
quency (Figs. 2, 3 in ref. II Oil , this volume), and V/r(r ms ), 
but these frequencies are fixed on human time-scales (the 
mass and angular momentum of the neutron star change 
only over millions of years). 

The resolution of the puzzle, suggested in ref. 
is non-linear resonance. Resonance occurs at character- 
istic frequencies, for non-linear systems these can vary 
00. 
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FIGURE 1. Radial v r , and vertical Vg epicyclic frequencies 
for a non-rotating black hole (spin a = 0) of mass M = lOM,-.,. 
Location of various possible resonances between the two is 
shown as a function of the dimensionless radius. 



EPICYCLIC RESONANCE 

Of course, one still needs to identify the resonant modes. 
An intriguing possibility is that the resonance occurs 
between accretion-disk motions occurring at the two 
epicyclic frequencies. An even more specific suggestion 
is that the resonance is parametric 1151 Hal , which can 
occur when the two frequencies are in a 2/n ratio, with 
n integer, fl3ll . Since the meridional epicyclic frequency 
is larger than the radial epicyclic frequency, V r /vg < 1, 
(Fig. 1), the lowest value of the integer at which para- 
metric resonance can occur is n = 3, assuming that it is 
driven by small-amplitude radial oscillations. This would 
give V r /Vg = 2/3 for the ratio of the two frequencies. A 
related g-mode calculation fl7ll yields 1 / v^2. 



SCO X-l 



The published data on Sco X-l J3 exhibits clustering of 
the kHz QPO frequency ratio near 2/3 1 19], suggesting 
that the twin QPOs may be related to parametric reso- 
nance between the epicyclic frequencies. However, the 
correlation line of the two frequencies has a different 
slope from 2/3. This can be explained if the system is 
slightly off -resonance 1 20] . 
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FIGURE 2. The correlation between the two kHz QPO fre- 
quencies in Sco X- 1 (after ref. 1 20]). The data (dots) is from van 
der Klis et al. 1 18]. The continuous line going through the main 
cluster of points is a theoretical calculation of the frequencies 
expected slightly off-resonance in a model of parametric reso- 
nance between epicyclic motions in an accretion disk. 



SIGNATURES OF NON-LINEARITY 

A rational, but not a unit, ratio of resonant frequencies, 
is expected in non-linear resonance, and the presence 
of QPO frequencies in a 2:3 ratio in microquasars l2lll 
is a strong argument in favor of this interpretation of 
QPOs. However, simple flute overtones can also give this 
ratio 12211 . see also 12311 - The presence of subharmonic 
frequencies is a hallmark of non-linear resonance 1 14]. 
Possible subharmonics have already been reported 1 24] 
in one of the microquasar QPOs, XTE J 1550-564, and 
their significance duly noted |16]. Strong evidence of 
non-linearity is afforded by the subharmonic (in relation 
to the neutron-star spin) separation of QPO peaks of 
about 200 Hz in the 401 Hz accreting pulsar 1 3]. 



FORCED RESONANCE 

A spinning neutron star at the center of an accretion disk 
is a source of possible perturbations for the latter. These 
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FIGURE 3. Peak power in the induced vertical oscillations 
of a slender torus perturbed radially at frequency v s , as a 
function of the ratio Vj/Av, where Av is the difference between 
the vertical and radial epicylic frequencies at the center of the 
torus. Three peaks are clearly visible, at 1:1, 3:2 and 2: 1 ratios. 



could be produced either through coupling via the pulsar 
magnetic field, or some disturbance on the surface of the 
star. In each case, the perturbation would be periodic, re- 
peating at the pulsar spin frequency, v v . Given this fact, 
one can ask which modes will be excited in the disk. Dy- 
namical (SPH) calculations of the oscillations of slender 
torii, slightly perturbed away from hydrostatic equilib- 
rium, show that the response to such perturbations occurs 
primarily at both, the local radial, and vertical, epicylic 
frequencies, even if the applied perturbation is purely ra- 
dial (or vertical) 12511 . This demonstrates that a certain 
amount of coupling, purely due to pressure, is present. 
We find 1 26], that the magnitude of the response is depen- 
dent on the relation between the spin period and the dif- 
ference between the epicyclic frequencies, Av. Peak re- 
sponse in vertical oscillations occurs when the difference 
between the two epicyclic frequencies is equal to one- 
half the perturbing frequency (Figure [3}- This, V s /2, is 
the frequency separation observed in the accreting 2.5 ms 
pulsar SAX J1808.4-3658 d. Other excitations, notably 
when the difference of epicyclic frequencies is equal to 
the pulsar spin frequency, are also possible. This would 
explain why in many neutron-star kHz QPO sources, the 
two QPO frequencies are separated either by one or one- 
half the inferred spin frequency 1 27|] . 

In conclusion, we concur that quasi-periodic modula- 
tions of the X-ray flux in low-mass X-ray binaries reflect 
oscillations of the accretion disk 128112911 . However, non- 
linear effects must be taken into account. 
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